Finite element analysis of fretting crack propagation 
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In this work, the finite elements method (FEM) is used to analyse the growth of fretting cracks. 
FEM can be favourably used to extract the stress intensity factors in mixed mode, a typical situation 
for cracks growing in the vicinity of a fretting contact. The present study is limited to straight cracks 
which is a simple system chosen to develop and validate the FEM analysis. The FEM model is tested 
and validated against popular weight functions for straight cracks perpendicular to the surface. The 
model is then used to study fretting crack growth and understand the effect of key parameters such 
as the crack angle and the friction between crack faces. Predictions achieved by this analysis match 
the essential features of former experimental fretting results, in particular the average crack arrest 
length can be predicted accurately. 

PACS numbers: 60.20Mk, 62.20Qp 
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I. INTRODUCTION 

Fretting damage is still a major issue for a large num- 
ber of multi-parts industrial applications subjected to vi- 
brations Many palliatives have been developed over 
the years such as modifying the geometry, coating sur- 
faces to reduce friction or using shot penning in some 
cases Q. Very small oscillatory displacement of surfaces 
in contact induce partial slip conditions which can lead 
to rapid crack initiation and imped the life of assem- 
bled structures. Well known examples can be blade/disk 
TA6V contact in jet engines or Al alloys riveted lap joints 
in aerospace structures. 

Various models have been developed over the years 
[3, El to predict the growth of fretting fatigue cracks but 
they usually fail to predict the early stages of the crack 
growth (typically before 1 mm). This stage is termed 
short crack behaviour and usually exhibits strong devi- 
ations with respect to the long crack domain which can 
be attributed to a combination of mixed mode growth 
under the influence of the highly multi-axial stress field, 
plasticity and roughness induced crack closure and crys- 
tallographic effects 0-0] ■ 

Recent developments in short fatigue crack models 
tend to take into account more and more of those in- 
gredients, but fretting crack growth is very often still 
described at the macroscopic level. Notable exception 
are recent studies which take into account the grain mi- 
crostructure in the contact region but they are limited 
to the initiation of fretting cracks [H, Q and a study of 
short fretting crack propagation using dislocation distri- 
butions • This work is a first step to study how FEM 
analysis can be favourably used to access local material 
behaviour at the crack tip to predict the overall fretting 
crack growth. We focus on the simple case of straight 
cracks initiated at the contact border; those cracks typi- 
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cally evolve in a multiaxial stress state which will strongly 
depend on the angle of propagation. Previous study of 
fretting crack propagation usually neglect mode II and 
use the assumption of Elbcr (AK = K max ) to determine 
the stress intensity factor range [H|,[llj]. While this might 
be justified for a fatigue crack under uniaxial load with a 
negative stress ratio, the case of fretting loading is much 
more complicated and there is no obvious reason for this 
assumptions to stand. In particular, the friction between 
the crack faces may influence mode II stress intensity 
levels due depending on the crack angle [13j . 

First the FE fretting model with straight cracks is de- 
scribed and two popular weight functions are used to 
test and validate the stress intensity levels in the case of 
cracks growing perpendicularly to the surface. FE com- 
putations are used to predict the mixed mode stress in- 
tensity levels during fretting cycles (thus accessing both 
K m ax and K m i n ) and to determine the influence of key 
parameters such as the crack angle and the effect of fric- 
tion between crack faces. Stress intensity ranges are even- 
tually combined to the Paris law to predict crack growth 
and arrest levels. Those predictions are compared to ex- 
perimental crack lengths observed on interupted tests in 
an 2xxx aluminium alloy. 



II. NUMERICAL METHODS 



A. Fretting model 

The simulations ultimately aim at predicting real fret- 
ting test and thus have been designed to closely match 
the experimental geometry, materials and loading condi- 
tions. 

The cylinder /plane fretting wear test shown in Fig.Uk) 
is represented in a 2D plane strain FE computation using 
the finite element code Z-SeT/ZcBuLoN [15| . Far from 
the contact zone, free meshing with plane strain 3-nodc 
triangular full integration elements is used and the mesh 
is gradually refined towards the contact region. To com- 
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FIG. 1: a) Simplified schematics of a fretting wear test in 
cylinder plane configuration, b) Measured fretting crack depth 
as a function of the number of cycles in a 2xxx series alu- 
minium alloy, P = 440 N/mm and Q = 240 N/mm; the three 
domains labelled 1,2 and 3 correspond respectivel y to crack 
initiation, fast propagation and crack arrest, after [14 ] . 



bine a high precision with reasonable computation times, 
the mesh in the contact zone is paved with plain strain 
rectangular 4-node full integration elements with a size of 
20 jim. With the investigated conditions, this provides a 
contact element over contact area ratio of 80 which is suf- 
ficient to accurately describe the multiaxial field imposed 
to the sample by the contact conditions (For a specimen 
free of cracks, the stresses computed with this geometry 
fits the analytical solution within 2%). 

The loading is applied in two steps: first, a normal dis- 
placement U2 is imposed on the top of the half cylinder to 
cause a normal reaction of 410 N/mm (corresponding to 
a half contact width of a — 0.8 mm and a Hertz pressure 
of 325 MPa). The plate is kept fixed by locking its bot- 
tom and side nodes. An horizontal cyclic displacement 
u\ is imposed to the bottom and side nodes of the plate 
to cause a tangential reaction of 240 N/mm. 

For the present study the material is described by 
isotropic elasticity; elastic constants E = 72 GPa and 
v = 0.3 are used to simulate a 2xxx aluminium alloy. 
Contact between the two parts is modelled by a classi- 
cal impactor/target technique and solved with a flexi- 
bility method implemented in Z-SeT/ZcBuLoN. Techni- 
cally, the contact reactions are at first computed in a local 
contact algorithm and then added to the global problem. 
The friction behaviour is introduced by the Coulomb's 
law with a friction coefficient of 1 . 1 , which represents the 
friction behaviour of a dry aluminium/aluminium con- 
tact [Hj]. 



contact area to receive the cracks. While both crack faces 
are initially superimposed, they are can open and close 
during the computation and the contact conditions be- 
tween them have to be defined. The crack tip is meshed 
with a fine, regular region to ensure a good description of 
the stress singularity; this region is used to compute the 
energy release rate G. All these aspects are schematically 
summarised in Fig. 2. 
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FIG. 2: Cracks are explicitly introduced within the finite ele- 
ments model with a refined mesh at the crack tip to compute 
the energy release rate G. 

Once cracks have been introduced within the mesh, 
finite element computation can be carried out to deter- 
mine the stress distributions in the neighbourhood of the 
crack tip during the fretting cycle. Due to the cyclic reci- 
procating motion of the counter body, both crack will al- 
ternatively open and close. As an example, Fig. 3 shows 
the mode I stress intensity factor for the left crack during 
one fretting cycle. As expected the left crack opens when 
the counterbody is moved to the right. By computing a 
complete fretting cycle, the stress intensity factor ranges 
AK] and AKjj can be determined. 

This analysis is used iteratively to compute the stress 
intensity levels during the growth of fretting cracks. Each 
different stage uses a different mesh accounting for the 
appropriate crack length and angle. The complete anal- 
ysis is as follows: 



B. Finite element analysis of fretting cracks 

In this approach, the cracks are explicitly introduced 
within the FE mesh. A in-house program is used to mesh 
both cracks independently, using as input their geome- 
try (here initial position of the crack, length and angle). 
Two rectangular spaces are reserved on both side of the 



1. For the first stage, the crack initiation position Xq 
and initial angle 9q have to be defined. This can 
be done using a fatigue multi-axial criterion (such 
as Dang Van or SWT) as discussed thouroughly by 
many authors (see [17| for a review). The present 
work focuses on propagation and those values are 
chosen somewhat arbitrarily, xq is chosen equal 
to 0.96a and 8q will vary from zero (crack is per- 
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FIG. 3: Example of left crack stress intensity factor extracted 
form the FEM calculation; here the crack is straight, 400 um 
long and inclined towards the center of the contact (9 = 15°). 



pendicular to the surface) to 45° (crack is inclined 
towards the center of the contact). 

2. The cracks are introduced within the mesh as ex- 
plained earlier and a complete fretting cycle is com- 
puted. 

3. The energy release rate G at the crack tip is deter- 
mined by post-processing the finite element results 
with a perturbation method [la ]. 

4. Local displacements fields u n and Ut respectively 
normal and tangent to the crack tip direction are 
computed and the ratio q — u t / u n is evaluated over 
the 4 pair of nodes behind the tip (this correspond 
to a 20 um distance). 

5. Both mode I and mode II stress intensity factors 
are deduced from 



and 



K u jKi 



with E' = E/(l — v 2 ) in plane strain. 

6. The crack geometry is updated with a new segment 
(li, 9i), 9i being kept constant and lo being equal to 
25 um. 



C. Weight functions for straight cracks 
perpendicular to the surface 

It is worth noting that thanks to R. D. Mindlin, the 
analytical solution of the partial slip regime is known 
precisely [191 ] . The analysis was derived for sphere/sphere 



contact but can be used also for a 2D cylinder/plane 
contact and provide the stress state of any material point 
M(x,y) of the plane during the complete fretting cycle 
(P held constant, Q(t) varying from +Q to —Q). 

For very small displacement amplitudes, slip occurs 
only on the outer edges of the contact, the central zone 
|x| < c remaining in stick condition. During the fretting 
cycle, the tangential force goes from +Q to — Q and the 
stick/slip limit d evolves from c to a. The partial slip 
contact can be described during the fretting cycle, by 
superimosing the elastic solution from the normal load 
only a n and three slipping contact elastic solutions of: 
(i) a slipping contact between —a and +a with a slip 
amplitude of +6 (ii) a slipping contact between — c and 
+c with a slip amplitude of +<5 (hi) a slipping contact 
between — d and +d with a slip amplitude of —28. 

Using this description, the stress state can be derived 
for any material point within the plane for a given fretting 
loading (P,Q). One can then describe the stress evolution 
along a virtual crack path during the fretting cycle [2fj| . 
In particular, at the end of the fretting cycle d — c the 
stress tensor at the position normalized by the conact 
width (x/a,y/a) is: 
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norm, load tang, forward load tang, reverse load 

Stress from Eq. ([I]) can be favourably used as input for 
a weight function of a straight crack to predict the stress 
intensity factor of fretting cracks. 

Many weight functions have been develo ped over the 
years since the seminal work of Bueckner [21|. In this 
work, two weight functions for a single edge crack in a 
plate will be used to validate the finite element approach 
in the case of a straight crack perpendicular to the sur- 
face. The first one is Bueckner weight function and the 
second one is an approximated weight function from Wu 
and Carlsson which account for the finite width of the 
plate [22{ . The stress intensity factor is expressed as 



K = fVnaW 
a is the crack length, W the specimen thickness and 



f = \W I a{x)m{a,x)dx 



(2) 



with ct{x) being the stress normal to the crack face along 
the crack path evaluated in the non cracked specimen 
from Eq. (J]) and 



i-3/2 



i=l 



where the Pi (a) functions can be found in [2z 
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III. RESULTS 

The model presented in section [TT] has been used to 
study the influence of the crack path amid the complex 
fretting stress field. First, straight cracks perpendicu- 
lar to the surface are investigated and weight functions 
are used to compare with values extracted from finite ele- 
ments calculations. Then the influence of the crack angle 
with respect to the surface on the stress intensity factors 
range evolution as the crack propagates is regarded; the 
effect of friction on crack faces is also studied. Eventu- 
ally life predictions are derived from the stress intensity 
factor ranges and compared to experimental values. 



FE prediction very well. Nevertheless this approach vali- 
dates the FE computations and the stress intensity levels 
predicted by this method. Regardless of the tangential 
force value, Kj first increases with crack length, reaches 
a maximum around 100 pm and then drops progressively. 
The drop in crack driving force was expected since with 
pure fretting conditions there is no bulk load and the 
stress level decreases when moving away from the con- 
tact area. The effect of the tangential force magnitude 
is also clearly highlighted by Fig. [5J increasing Q lead to 
an increase in maximum crack driving force. The mag- 
nitude of Q also seems to influence the drop of Ki, and 
one can qualitatively predict a longer crack with a higher 
Q which is in agreement with experimental observations. 



A. Straight crack perpendicular to the surface 

Cracks are introduced in the mesh at xq = ±0.77. Both 
stress intensity factor ranges are extracted at each time 
increment during the fretting cycle. After each calcu- 
lation (starting with a 25 pm crack), the crack length is 
increased by 25 pm. Selected results are shown on Fig. [4] 
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FIG. 4: Isovalues of von Mises stress for 3 different crack 
lengths: 50, 150 and 450 um; straight crack perpendicular to 
the surface. 
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FIG. 5: Comparison of mode I stress intensity factors for a 
straight fretting crack perpendicular to the surface computed 
by means of weight functions and finite elements. 




Computations have been carried out for two different 
values of the tangential force (Ql — 150 N/mm and 
Q2 = 240 N/mm). The maximum stress intensity fac- 
tor for mode I has been extracted as explained in sec- 
tio rill Bl and are compared to the values predicted by the 
use of weight functions (see section III C[) . For the case 
of cracks perpendicular to the surface, normal loading of 
the crack face is a xx (y) described by which is obtained 
analytically from Mindlin analysis and used in Eq. @ 
to derive the evolution of Kj as a function of the crack 
length. For small crack lengths, the values computed by 
FE are in very good agreement with those calculated with 
weight functions as for longer crack lengths, the value 
computed by FE is slightly lower. Presumably this can 
be attributed to the finite length of the sample which is 
not taken into account in the weight functions used here; 
a similar effect can been observed in Single Edeg Notched 
Tensile (SENT) specimens under uniaxial loading. In- 
deed, computing Kj by Bueckner's weigth function de- 
viates from the FE values when the crack size increases. 
For the SENT case, a solution accounting for the finite 
size of the specimen [23[ can be used which match the 



B. Slanted crack 

Experimental observations of fretting cracks by met- 
allography techniques often show that cracks can be in- 
clined towards the center of the contact. The angle is 
expected to have an influence on the mode I and mode 
II levels due to the multiaxial stress state in the subsur- 
face contact region. FE computations series have been 
carried out with several values of the crack angle ranging 
from 15° to 45° (0° being perpendicular to the specimen 
surface). Selected results are shown on Fig. [5] for the 
cases of 30°. 

The stress intensity factors have been computed as ex- 
plained in section III Bl for a variety of crack angles and 
crack lengths. A comparison for the cases of and 30° is 
given in Fig. [7] for different crack lengths. 

Fig. [7] highlights the importance of considering the 
whole fretting cycle. For Kj levels, the maximum is al- 
ways reached at the end of the fretting cycle, where the 
crack is open. For cracks perpendicular to the surface, 
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FIG. 6: Isovalues of von Mises stress for 3 different crack 
lengths: 50, 150 and 450 um; straight crack inclined with a 
30° angle. 
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FIG. 7: Energy release rate G and stress intensity factors 
Ki and Ku evaluated during the fretting cycle: left column 
shows results for a straight crack perpendicular to the surface 
whereas right column shows results for an inclined crack (8 — 
30°); period of the fretting cycle is 8 seconds. 



intensity factor range AKj as a function of crack depth 
for all investigated crack angles (0, 15, 30 and 45°). 
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FIG. 8: Evolution of the stress intensity factor ranges AKj 
and AKn as a function of crack depth for all investigated 
crack angles. 



Effect of friction on crack faces 



Ki min remains close to zero and the decrease in terms 
of stress intensity factor range is due to the decrease of 
Ki max . For slanted cracks, the drop of AKi as a func- 
tion of crack depth is more pronounced due to the crack 
angle. The more the crack is oriented towards the cen- 
ter of the contact, the faster the crack driving force will 
drop. Considering Ku levels, the maximum occurs at a 
different time during the fretting cycle for crack perpen- 
dicular to the surface and slanted cracks. It is interesting 
to see that for slanted cracks, Kjj max values will not de- 
crease to zero as for Ki mam , but increase. However, Ki min 
also increases, resulting in a decrease of AKn levels as 
a function of crack depth are shown on Fig. [5] We shall 
demonstrate in the next section that this can be predom- 
inantly attributed to the friction of the crack faces. 
Fig. [5] shows a summary of the evolution of the stress 



To check the influence of friction applied between the 
crack faces, computations for the crack inclined with a 
30°angle has been carried out a second time enforcing 
zero friction between the crack faces. Results of these 
calculations are depicted in Fig. [9] This does not change 
the maximum values of the stress intensity factors, since 
those are obtained at the end of the fretting cycle when 
the crack is open. However, this clearly shows that the 
minimum value of Ku is significantly lower leading to 
higher AKjj levels. This is due to the fact that slanted 
cracks evolve in mixed mode and friction will limit crack 
tip shear displacement. If friction of crack lips is en- 
forced, stick condition will be reached at some point pre- 
venting the shear displacement of crack tip nodes and 
leading to a Kjj min value close to Kjj max as seen on 
Fig. [7] for a crack larger than 350 um. This result will 
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strongly impact life predictions (maximum crack depth 
as well as number of cycles to reach a certain depth) when 
considering mixed mode crack growth. 




50 100 150 200 250 300 350 400 450 500 
crack depth ([jm) 



D. Life Predictions 

Using the previous calculations, life predictions based 
on the computed stress intensity factor range -or more 
exactly crack depth prediction as a function of the num- 
ber of cycle, since specimens do not fail in pure fretting 
conditions- can be derived from the Paris Law: 



N = Ni + 



1 



io CAR 



-da 



eff 



where K P 



y/ + Kj T , the coefficient C and the ex- 
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ponent m, are given for da/dN expressed in mm/cycle 
and K in MPa-y/m. This equation has been originally 
developed for mode I fatigue crack growth but can also 
reasonably well describe mixed mode crack growth as 
shown in 26]. The two constants are equal to C = 
6.53 x 10~ 8 mm/cycle, m = 3.387 for the considered 2xxx 
aluminium alloy. For the studied fretting conditions, Ni 
has been experimentally determined to be equal to 50. 10 3 
cycles [lj|- 
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FIG. 9: Effect of the friction between crack faces on mode II 
stress intensity levels. 



Those results highlight the key role of the friction be- 
tween crack faces. This effect is particularly expected in 
the absence of bulk fatigue load (in pure fretting con- 
ditions) where the mode I stress intensity factors will 
progressively decrease to zero and mode II will become 
predominant. The crack progression will then strongly 
be controlled by the friction between the crack faces. A 
high friction will prevent shear tip displacement and stop 
the crack. High observed friction in experiments could 
also be related to roughness induced crack closure, since 
real fretting cracks may exhibit either strong deviations 
due to a crystallographic crack path, either small rough- 
ness as observed by post mortem examination or X-ray 
tomography slices [2J] . In both case, RICC levels are ex- 
pected to be quite high which would lead to lower AK 
levels and a lower crack arrest depth '2&\ . 



FIG. 10: Life predictions based on the stress intensity factors 
determined by FEM for different crack angles; crack depth 
measured experimentally have been added for comparison. 

Fig. QUI shows the results of the fretting life predictions 
for the same loading conditions as used in the experi- 
ments [3| • One can readily observe that the crack angle 
has a strong influence on the propagation life as sug- 
gested by the AKi and AKjj values already presented. 
Increasing the angle of the crack towards the center of 
the contact will have two effect: (i) the crack growth will 
be slower and (ii) the crack arrest depth will be shorter. 
For clarity, predictions for the 30° angle without fric- 
tion between the crack faces are not presented in Fig. [10] 
but the curve is essentially the same with a longer crack 
arrest depth (300 um versus 200 um in this particular 
case). This is clearly due to the higher AKu value which 
rapidly becomes predominant for inclined cracks (AKj 
reaches zero). 

Regarding the comparison with the experimental re- 
sults, one can observe that the 30° angle predictions 
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closely match the initial experimental behaviour, al- 
though the average crack arrest behaviour is best de- 
scribed by the 15° angle results. This is to be compared 
with the 20 to 25° average propagation angle measured 
on cross section and by tomography. Of course in real 
experiments, the crack angle is not constant and further 
work is needed (using a bifurcation criterion) to extend 
those prediction with a varying crack angle as a func- 
tion of the crack depth. Criteria for the determination 
of angles of crack growth are now well established un- 
der proportional loading (27l - [29| . However, fretting con- 
ditions induces non-proportional loading at crack tips 
which arc characterized by a Kj / Kn ratio varying dur- 
ing the cycle. Only a few theorical studies are devoted 
to non-proportional conditions and some of them [30l l3l| 
showed the inapplicability of certain classical criteria to 
describe crack propagation. For instance, Hourlier and 
Pineau applied a criteria under non-proportional condi- 
tions where the crack path corresponds to the direction 
where the mode I fatigue growth rate is maximum [32| . 
Bower suggested that the crack would tend to propagate 
in the direction of the maximum value of Aa max }33| . 
Lamacq showed that for fretting fatigue conditions, the 
Kii = approach worked well to predict the crack an- 
gle (26[. Yet, there is no general agreement on which 
criterion is to work in the pure fretting case. 

IV. CONCLUSION 

This paper presented a numerical framework suitable 
to study crack propagation under pure fretting condi- 
tions. The FEM simulations have been coupled with 



LEFM to successfully describe some experimental char- 
acteristic features of fretting cracks like the influence of 
the crack angle and the crack arrest depth. Given the 
highly multi-axial stress field developing under the con- 
tact, the path taken by the crack appears critical with 
respect to its propagation rate. Regarding mode I, the 
maximum stress intensity value decreases as the crack 
grows and it was shown that the decrease is stronger 
when the crack is inclined towards the center of the con- 
tact zone. For mode II, although the stress intensity level 
Kii remains high, its range AKn progressively decrease 
to zero if the friction of crack faces is taken into account, 
this ultimately being the cause of the crack arrest. If 
crack propagation is assume to occur without friction, 
predicted mode II stress intensity factor ranges are much 
higher which leads to longer fretting cracks. 

Two directions appear particularly important to ex- 
tend this work. First, the effect of crack bifurcation 
should be investigated. Many criteria are available from 
the literature to try to predict the path of a crack in 
a multiaxial stress field but the non-proportional nature 
of fretting loads should be accounted for. Second, the 
present model cannot yet account for microstructural ef- 
fect on the crack path as shown by a number of stud- 
ies. Cyclic plastic deformation at the crack tip should 
be investigated using a suitable material constitutive be- 
haviour such as crystal plasticity. 

Acknowledgments 

The authors are grateful to Professor G. Cailletaud for 
fruitful discussion and comments. 



[1] R. B. Waterhouse, International Metals Reviews 37, 77 
(1992). 

[2] S. Fouvry, V. Fridrici, C. Langlade, P. Kapsa, and L. Vin- 
cent, Tribology International 39, 1005 (2006), ISSN 
0301-679X, the Fourth International Symposium on Fret- 
ting Fatigue. 

[3] M. Szolwinski and T. N. Farris, Wear 221, 24 (1998). 

[4] J. Madge, S. Leen, and P. Shipway, International Journal 
of Fatigue 30, 1509 (2008), ISSN 0142-1123. 

[5] D. L. McDowell, Engineering Fracture Mechanics vol. 
56, 357 (1996). 

[6] V. Doquet, Fatigue and Fracture of Engineering Materi- 
als and Structures vol. 22, 215 (1999). 

[7] W. Ludwig, J.-Y. Buffiere, S. Savelli, and P. Cloetens, 
Acta Materialia 51, 585 (2003). 

[8] R. Neu, Key Engineering Materials 378-379, 147 (2008). 

[9] T. Dick, S. Basseville, and G. Cailletaud, Computational 
Materials Science 43, 36 (2008), ISSN 0927-0256, pro- 
ceedings of the 16th International Workshop on Compu- 
tational Mechanics of Materials - IWCMM-16. 
[10] S. Fouvry, D. Nowell, K. Kubiak, and D. Hills, Engineer- 
ing Fracture Mechanics 75, 1605 (2008), ISSN 0013- 
7944. 



[11] H. Fadag, S. Mall, and V. Jain, Engineering Fracture 
Mechanics 75, 1384 (2008). 

[12] E. Giner, N. Sukumar, F. Denia, and F. Fuenmayor, In- 
ternational Journal of Solids and Structures 45, 5675 
(2008). 

[13] R. Ribeaucourt, M.-C. Baietto-Dubourg, and 
A. Gravouil, Computer Methods in Applied Me- 
chanics and Engineering 196, 3230 (2007), ISSN 
0045-7825. 

[14] S. Munoz, H. Proudhon, J. Dominguez, and S. Fouvry, 

International Journal of Fatigue 12, 1769 (2006). 
[15] J. Besson and R. Foerch, Revue Europeenne des 

Elements Finis 7, 535 (1998). 
[16] H. Proudhon, S. Fouvry, and J.-Y. Buffiere, Inter- 

national Journal of fatigue 27, 56 9 (2005), URL 

|http : I It r . arxiv . org/abs/0710 ■ 3715] 
[17] S. Fouvry, K. Elleuch, and G. Simeon, Journal of strain 

analysis 6, 549 (2002). 
[18] D. M. Parks, International Journal of Fracture 10, 487 

(1974). 

[19] R. D. Mindlin and H. Deresiewicz, Journal of Applied 

Mechanics vol. 20, 327 (1953). 
[20] D. A. Hills and D. Nowell, Mechanics of Fretting Fa- 



8 



tigue, 0792328663 (Kluwer Academic Publishers, Dor- 
drecht, Boston, 1994), ISBN 792328663. 

[21] H. F. Bueckner, Zeitschrift fr Angewandte Mathematik 
und Mechanik 50, 529 (1970). 

[22] X. Wu and J. Carlsson, Weight functions and stress in- 
tensity factor solutions (Pergamon Press, 1991), ISBN 
0080417027, 9780080417028. 

[23] J. Ahmad, V. Papaspyropoulos, and A. T. Hopper, En- 
gineering Fracture Mechanics 38, 283 (1991). 

[24] H. Proudhon, J.-Y. Bufnere, and S. Fouvry, Engineering 
Fracture Mechanics 74, 782 (2007). 

[25] M. R. Parry, S. Syngellakis, and I. Sinclair, Materials 
Science and Engineering A 291, 224 (2000). 

[26] V. Lamacq, M. C. Dubourg, and L. Vincent, Journal of 



Tribology 118, 711 (1996). 
[27] F. Erdogan and G. Sih, J. Basic Engng. 85, 519 (1963). 
[28] G. Sih, International Journal of Fracture 10, 305 (1974). 
[29] J. Tirosh, Engineering Fracture Mechanics 9, 607 (1977), 

ISSN 0013-7944. 
[30] P. E. Bold, M. W. Brown, and R. J. Allen, Fatigue and 

Fracture of Engineering Materials and structures 15, 965 

(1992). 

[31] E. W. Smith and K. J. Pascoe, Fatigue and Fracture of 
Engineering Materials and structures 6, 201 (1983). 

[32] F. Hourlier and A. Pineau, Fatigue of engineering mate- 
rials and structures 5, 287 (1982). 

[33] A. Bower, ASME Journal of tribology 110, 704 (1988). 



